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ABSTRACT 

Considerable interest is now focused on the detection of terrestrial mass planets around M dwarfs, and radial 
velocity surveys with high-resolution spectrographs in the near infrared are expected to be able to discover 
such planets. We explore the possibility of using commercially available molecular absorption gas cells as a 
wavelength reference standard for high-resolution fiber-fed spectrographs in the near-infrared. We consider the 
relative merits and disadvantages of using such cells compared to Thorium- Argon emission lamps and conclude 
that in the astronomical H band they are a viable method of simultaneous calibration, yielding an acceptable 
wavelength calibration error for most applications. Four well-characterized and commercially available stan- 
dard gas cells of H'^C'^^N, '^C2H2, '^CO, and '^CO can together span over 120nm of the H band, making them 
suitable for use in astronomical spectrographs. The use of isotopologues of these molecules can increase line 
densities and wavelength coverage, extending their application to different wavelength regions. 

Subject headings: techniques: radial velocities — techniques: spectroscopic - 
— stars: kinematics — methods: data analysis 



- instrumentation: spectrographs 



1. INTRODUCTION 

With the discovery of more than 300 extrasolar planets, 
considerable interest is now focused on finding and char- 
acterizing terrestrial mass planets in habitable zones around 
their host stars. Such planets are extremely difficult to detect 
around F, G and K stars, requiring either extremely high ra- 
dial velocity precision (< 1 m s"') or space based photometry 
to detect a transit. Radial velocity and transit efforts are now 
beginning to focus on M dwarfs, the most numerous stars, 
where the lower luminosity shifts the habitable zone much 
closer to the star Terrestrial mass planets in such regions are 
detectable with the current radial velocity precision obtained 
with high resolut ion echelle spectrographs (iButler et alJ 19961 : 
iPepe et aT]|2004l) . A number of the brighter M dwarfs are be- 
ing surveyed in the optical with existing high precision radial 
velocity instruments. Neptune and super-Earth mass planets 
have been discovered around a number of such objects. Ex- 
amples are GJ436 (Buder et al. 2004), GL581 (Udry et al. 
2007), and GL176 (Forveille et al. 2008) and they suggest 
that such planets may be rather common around M stars. The 
habitabilit y of terrestr i al pla nets around M stars has been ex- 
plored bv Tarter et al.1 (l2007h and lScalo et akl (l2007 b. 

However, most of these stars are intrinsically faint in the 
optical, emitting most of their flux in the 1-1.8 fim wave- 
length region, the near infrared (NIR) J (L 1-1.4 fim) and H 
(1.45-1.8 /im) bands. Stellar activity, coupled with the rel- 
ative faintness, can make detections in the optical difficult. 
For example, Endl et al. (2008) claim the detection of a 24 
Earth-mass planet around Gl 176 using 28 radial velocities 
obtained with the High Resolution Spectrograph (HRS) on the 
Hobby-Eberley Telescope, while Forveille et al. (2008), with 
58 higher precision measurements obtained with the HARPS 
(High Accuracy Radial velocity Planet Searcher) instrument, 
find evidence for stellar activity as well as a different orbital 
period and a much smaller minimum mass. Similar observa- 
tion in the NIR may require less telescope time and are less 
prone to systematics because the activity induced radial ve- 
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locity jitter is expected to be smaller in the NIR than in the 
optical. As an example, Setiawan et al. (2008) announced the 
discovery of a planet around the young active star TW Hya, 
but Huelamo et al. (2008) see no such variability on preci- 
sion radial velocity data obtained with the CRJRES (Cryo- 
genic Infrared Echelle Spectrograph) instrument in the infra- 
red (Seifhart & Kaufl 2008) and they attribute the velocity 
variability in the optical to star spots. Rucinski et al. (2008) 
also observe photometric periodicity at the proposed planet's 
period in one season of observations with the MOST (Mi- 
crovariability & Oscillations of Stars) satellite, and the varia- 
tions are absent in another season- suggesting that the cause 
of the radial velocity is activity-induced. Prato et al. (2008) 
have discussed radial velocity observations of young stars in 
the optical and NIR and conclude that observations in the 
near-infrared are essential to discriminate between activity 
and the presence of planets. NIR spectroscopy can explore 
new regimes of planet formation as well as complement ex- 
isting observational programs in the optical that target young 
or active stars. 

In this article, we explore the calibration challenges faced 
by a high-resolution spectrograph operating in the NIR and 
suggest that easily available absorption gas cells can be used 
for wavelength calibration purposes for any such instruments 
that are currently being designed. Our discussion focuses on 
using a simultaneous wavelength calibrator along a separate 
optical fiber (Baranne et al. 1996) to track and calibrate out 
instrument drift, similar to the Th-Ar calibration technique in 
the optical. We briefly highlight important aspects on such in- 
struments in §2. In §3 and §4 we discuss currently used emis- 
sion wavelength sources and their relative merits and disad- 
vantages. In §5 we explore the use of commercially available 
absorption cells for precision wavelength calibration. 

2. FIBER-FED HIGH RESOLUTION SPECTROGRAPHS IN THE NIR 

The use of optical fibers allows a spectrograph to be placed 
in a stable environmentally controlled enclosure. The intrin- 
sic radial scrambling properti es of an optical fiber can , with 
the use of a double scrambler (iHunter & Ramseylll992h . lead 
to a very stable illumination profile on the spectrograph slit. 
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Commercially available optical fibers provide high transmis- 
sion in the J and H NIR bands. Fiber-feeds allow a second 
calibration fiber to be used to simultaneously track the instru- 
ment drift during an object exposure. If the science and cal- 
ibration fiber images are relatively close on the focal plane, 
then the measured drift of the calibration spectrum is a very 
good estimate of the drift estimated on the science spectrum. 
The major requirements on the calibration source are that it be 
stable enough to achieve the required precision, have a num- 
ber of lines with known wavelengths, and be bright enough 
that sufficient S/N is obtained for it to not compromise the 
possible radial velocity accuracy on bright sources. Design 
studies for such instruments in the NIR have been conducted, 
including one for a proposed Precisio n Radial Velocity Spec- 
trograph ( iRavner & PRVS TeamI 120071) for the Gemini tele- 
scope. ^On_the_^ber-fedbench top PRVS Pathfinder instru- 
ment, iRamsev et al.l (l2008l) have demonstrated a short term 
radial velocity precision of 7-10 m s"' in the NIR using in- 
tegrated sunlight with Thorium- Argon as a wavelength ref- 
erence. The use of silicon immersion gratings is also being 
explored (eg. Ge at al. 2006, Jaffe et al. 2006) to facilitate 
compact spectrographs capable of providing spectral resolu- 
tions of R=50-100k in the NIR. 

3. THORIUM-ARGON EMISSION LAMPS: 

3.1. History & Advantages 

iKerber et alJ (120071) provide a summary of the development 
and use of Thorium-Argon (Th-Ar) hollow cathode emission 
lamps. The only naturally occurring Thorium (Th) isotope, 
^•'^Th, has zero nuclear spin, leading to sharp symmetric emis- 
sion lines, even at very high resolutions. The monatomic Ar- 
gon gas also has a number of bright lines. Th-Ar emission 
lines span the UV-NIR regions, making such a lamp a very 
useful and convenient source of wavelength calibration. Most 
wavelength calibration applications using this lamp are tied 
to the Palmer & Engleman ( 1983) (PE83 hereafter) measure- 
ments of the Thorium lines in the 280-1 lOOnm region. Their 
quoted measurement precision, ^ 0.001 cm"' - 0.005 cm"', 
corresponds to an intrinsic velocity uncertainty of 15-80 ms"' 
per emission line (at 550nm). For Argon (Ar) lines, the wave- 
length measurements typically used are adopted from Norli3 
(119731) . More recent measurements of Ar lines have been 
taken by Whaling et al. (1995, 2002), who extend measure- 
ments to the mid-IR. Argon is not a heavy element like Tho- 
rium and the central wavelengths of the Argon emission lines 
are susceptible to pressure shifts, making them sensitive to 
environmental conditions in the lamp. The lines may be un- 
stable by many tens of ms"' (iLovis & Pepell20()7l) . and for the 
highest possible precision (~ 1 - 3 ms"' ), Ar lines should gen- 
erally be avoided. At intermediate spectral resolutions, a fur- 
ther problem is line blends, because the intensity ratios of Th 
and A r lines is a function of the current supplied ( Kerbe r et al.l 
l2007h . Very stable spectrographs offer the opportunity of 
decreasing the internal errors associated with the original 
Thorium measurements and PE83s use of the Fourier Trans- 
form Spectrograph (FTS) at Kitt-Peak National Observatory 
(KPNO). Stable instruments allow one to acquire many Th-Ar 
spectra and co-add them, thereby increasing the S/N and re- 
ducing the photon noise uncertainties on the line-centers of 
known lines, as well as enabling wavelength estimation of 
lines not initially present in the or iginal PE 83 atlas. Such 
a technique has been applied by Lovis & Pepe (2007) with 
the vacuum enclosed HARPS high-resolution echelle spectro- 
graph, resulting in a line-list with improved wavelength pre- 



cision (though inheriting the zero-point of the original KPNO 
atlas). This improved line-list has certainly helped HARPS 
achieve radia l velocity precisio ns of <1 ms"' and discover 
super- Earths dMavor et al.ll2008l) . 

3.2. Calibration in the NIR 

In the NIR J and H bands Th lines are few and relatively 
faint, making it difficult to find accurate dispersion solutions 
in th e J and H band without substantial integration time. The 
iHinkle et al.l (120011) atlas contains FTS and grating spectro- 
graph wavelength measurements of 500 lines in the 1-2.5 fim 
range and Engleman, Hinkle & Wallace (2003) derive FTS 
wavelengths for a much larger sample using a cooled Th-Ar 
cell operated at a high current. Kerber et al. (2008) also mea- 
sured wavelength for Th-Ar lines for multiple cells with an 
FTS for use in calibration of the CRIRES instrument model. 
Although these atlases help to identify lines when they are de- 
tected, the Th lines in commercially-available lamps are still 
very low in intensity. The high contrast between the Th and 
Ar lines leads to the bright Ar lines dominating the spectra and 
these can be a source of significant scattered light. For exam- 
ple, the SOPHIE instrument (Perruchot et al. 2008) requires 
the use of a 700 nm shortpass filter to block scattered light 
from bright Ar lines in the infrared. For high precision ra- 
dial velocity studies, the situation is further exacerbated by the 
need to acquire a high enough signal level on sufficient lines 
to achieve a dispersion solution accurate enough to be able 
to measure radial velocities at a high precision. Another mi- 
nor concern is that simultaneous calibration with Th-Ar also 
makes scattered light subtraction of echelle frames difficult. 
The need to have an acceptable signal on most lines leads to 
the strongest lines saturating and bleeding into the adjacent 
stellar spectrum, making reduction more complex. The rela- 
tiv e lack of suffic i ent T h lines in the NIR has also been noted 
by Ramsey et aP (l2008l) in their laboratory tests with PRVS 
Pathfinder and these authors propose combining light from 
multiple lamps (most notable Uranium- Ar) with Th-Ar to in- 
crease the line density available for wavelength calibration. 

4. BROADBAND FREQUENCY COMBS 

The calibration advantages of using a laser frequency comb 
with high resolution spec trographs have be en discussed by 
iMurphv et akl (l2007h and iBraie et all (l2008h . Though these 
have the desired frequency coverage, the intrinsically small 
frequency mode spacing (250MHz- IGHz) makes their di- 
rect use unsuitable for astronomical spectrographs since the 
individual lines would blur together. The use of a Fabry 
Perot filter to sel ect only w ell-spaced emission lines has been 
demonstrated bv iLi et aU (2008) for optical wavelengths and 
in parts of the NIR H band (1530-1600nm)by Steinmetz et al. 
(2008), who achieved high-wavelength precision when using 
the comb for solar observations. Although there is little doubt 
that such femtosecond frequency combs offer the highest pos- 
sible precision, they are not yet easily available and are rela- 
tively expensive (though their cost is expected to come down 
in the future). While facility-class instruments on large tele- 
scopes will likely be able to invest in such calibrators, they are 
not as easy to obtain, maintain and operate as Th-Ar or other 
emission lamps. Many high resolution NIR spectrographs 
may also not need this level of sub 1 ms"' wavelength calibra- 
tion since their performance may be dominated by other sys- 
tematics. The atmospheric absorption and OH emission lines 
may also set the limit to the precision that can realistically 
be achieved in the NIR. Recent progress in generating THz 
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spacing frequency combs (Del'Haye et al. 2007) by using a 
micro-toroidal oscillator and a laser may offer some promise 
in making an inexpensive, readily available device for future 
NIR spectrographs. 

5. ABSORPTION GAS CELLS 

Absorption gas cells have a long history of usage for cali- 
bration. Molecular absorption caused by rotational and vibra- 
tional energy levels is we ll unde rstood, stable, and repeatable. 
The iCampbell & Walker! (11979 ) radial velocity survey used 
the sharp isolated absorption lines of HF gas in the red to cal- 
ibrate out instrument drifts. The use of molecula r Iodine (I2) 
gas as a simultaneous absorption calibrator by Butler et aH 
(119961) has lead to the ability to measure very precise veloci- 
ties, with long term precision approaching 1 ms"', and the sta- 
bility of such cells is now widely accepted. Absorption cells 
are not quite as popular as Th-Ar in the UV-optical region be- 
cause the lines tend to be localized in a smaller wavelength 
regimes (eg. 500-620nm for I2). In the case of I2, the lines are 
all blended even at spectral resolutions of 100k, requiring FTS 
spectra of that particular cell and PSF modeling to determine 
a wavelength scale when used with a spectrograph. Absorp- 
tion cells are in use in some astronomical instruments in the 
infra-red for wavelength calibration. The CRJRES instrument 
on the VLT uses a N2O gas cell as a calibrater for wavelengths 
longer than 2.2 /im (Kerber et al. 2008). Gas mixtures for an 
l2-like absorbtion gas cell for the planned N AHUAL spectro- 
graph have been discussed by Martin et al.l 12005) and for the 
planned GIANO spectrograph by D' Amato et al. (2008) who 
have developed a combined HCL-HBr -HI cell th a t has ~ 200 
lines, spanning J, H, & K bands. Ram sey et al.l ( l2008l) also 
mention plans to use HF and water vapour cells in their J-band 
pathfinder instrument to monitor velocity drifts. HF has deep 
lines, but they are few, and the gas is corrosive and difficult to 
work with. Such a cell also needs to be heated to more than 
70 °C to prevent polymerization of HF (Campbell & Walker 
Il979l) . resulting in line broadening. Water vapour also has 
many lines, but is also present in the atmosphere, making 
identification of isolated individual lines for wavelength cali- 
bration quite challenging. 

In the NIR H band (~ 1.45-1.8 /im) a number of absorp- 
tion calibrators with well spaced lines have been well char- 
acterized, primarily owing to the needs for ever higher wave- 
length division multiplexing in the telecommunication indus- 
try. Although no single cell covers the entire spectral region, 
we show that a significant fraction of the H band can be cov- 
ered using combinations of four commercially available cells. 
Our proposed use is not to pass the starlight through the cells, 
but to use the isolated absorption lines as wavelength mark- 
ers to determine the dispersion solution and track instrument 
drifts, thereby enabling the measurement of precise velocities 
and accelerations. 

5.1. NISTSRMS 

The United States National Institute of Standards & Tech- 
nology (NIST) has designated four standard reference mate- 
rials (SRMs hereafter) for use in wavelength calibration. To- 
gether these four materials span the C & L telecom wave- 
length bands and have ^ 200 isolated sharp lines in the 1510- 
1630 nm wavelength region of the H band. Table 1 lists these 
NIST SRM gas cells, their effective wavelength coverage and 
the number of lines with NIST-certified wavelengths. These 
gas cells are commercially available, and have an effective ab- 
sorption path length of 5-80 cm. The low pressure H'^C'^^N 



(Swann & Gilbert 2000) and '2C2H2 dSwann & GilberllllOOSh 
cells have been primarily designed for high resolution appli- 
cation and have narrow linewidths of 7-15 pm, essentially 
unresolved at spectral resolutions as high as R=50-70k at 
1550nm. The high pressure CO gas cells were originally de- 
signed for use with an inst rumental resolution of .05nm and 
have linewidths of 50pm ( iSwann & Gilbertll2002h . Decreas- 
ing the cell pressure can easily make the lines sharper, if nec- 
essary, and increasing cell length can increase absorption for 
unsaturated lines. These gas cells are available coupled to 
single-mode fibers with standardized FC/PC connectors for 
light input and output, allowing the CO gas cells to be used 
in a multi-pass configuration, achieving 80 cm of absorption 
length with four passes through a compact 20 cm cell. The 
use of the single mode fiber of input and output makes it easy 
to pass continuum light through multiple cells to obtain an im- 
print of the absorbtion lines over larger wavelength regimes. 
The use of the fibers is not necessary, but is a convenience. 

5.2. Astronomical Echelle Spectrographs 

Modern high resolution spectrographs generally use coarse 
ruled echelle gratings that have high blaze angles (R2-R4 or 
tan 6*5 = 2-4 is fairly typical). To estimate the achievable ac- 
curacy of wavelength calibration, we consider an R2 echelle 
spectrograph capable of a resolution of R=50k at 1600 nm. 
We assume a 4 pixel sampling the FWHM of the Gaussian 
resolution element. We assume a grating with a groove den- 
sity such that a wavelength of 1550nm corresponds to an or- 
der number of ^ 100. While such low groove densities are 
not available in commercially ruled echelles, they c an be cre- 
ated on silicon immersion gratings dGe et alj|2006l) enabling 
high dispersion as well as the ability to pack all orders from 
the H band into a single Ik by Ik detector. Figure 1 shows the 
expected absorption spectra from these gas cells with such an 
instrument. For HCN an d C2H2 we have u sed the high reso- 
lution spectra scans from lSwann & Gilbert! (HOOO. . 2005)^ and 
convolved them with a Gaussian instrument profile (R=50k at 
1600nm). For CO the full spectral scans are not available and 
we have simulated t he data using line cente rs and broaden- 
ing coefficients from ISwann & GilbertI (2002). We have not 
convolved the simulated CO data with the Gaussian instru- 
ment profile since the measurements were obtained at a lower 
resolution than R=50k. We would expect these lines to be 
slightly deeper when observed at the resolution we assume for 
the spectrograph. Figure|2]shows a schematic arrangement of 
the gas cells, coupled with optical fibers, that can be used to 
generate the absorption spectra seen in Figure 1 . 

The effective free spectral range of each order, ^ X/m, cor- 
responds to 15nm bandwidth at 1550nm. As can be seen from 
Figure 1, the HCN cell alone provides 20 sharp reference lines 
within ±8 nm of the order center, a region where th e blaze ef- 
ficiency is expected to be > 40% of peak dSchroe der 1987). 
In reality the echelle order spans a larger wavelength cover- 
age than ^ X/m on the CCD, making more lines accessible 
for wavelength calibration. The design of echelles also leads 
to some wavelength overlap between adjacent orders, allow- 
ing the same line to be used to calibrate more than one or- 
der Gratings operating at lower order numbers will have even 
higher line densities per order if the orders fit on the detector 
array. 

Typical echelle spectra span 7-12 pixels in the sUt direction, 
and we assume 7 pixels here. Assuming that the continuum 

^ Also available at www.nist.gov/srm 
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light source is bright enough (discussed in detail later), a typ- 
ical continuum S/N of 200 per pixel is a reasonable assump- 
tion, assuring also that one is operating in the linear regime in 
NIR detectors (which typically have quantum wells of 100k). 
The data reduction process will optimally combine all 7 pixels 
in the slit direction into a single pixel with a S/N of 200 \/7, or 
S/N ^ 500 per pixel, in the extracted one-dimensional spec- 
trum. Although the actual shape of the absorption line is a 
Voigt profile, we can approximate the line as a Gaussian for 
the purpose of estimating the velocity precision. Following 
the procedure outlined in iButler et aH (11996.) for calculating 
the photon-noise limited velocity precision yields an error of 
24 m s"' for a line with a depth (after convolution with the 
spectrograph instrument profile) of 20%. So the use of ^20 
lines reasonably well spaced across the echelle order can, in 
principle, enable wavelength calibration accuracy at the level 
of ^ 5 m s"' per echelle order This uncertainty is smaller 
than the photon noise error expected per order on most stars 
observed, and is therefore sufficient. Many HCN and C2H2 
lines are deeper than 30%, while the CO lines have lower 
depths of ~ 10-18% and have line widths broader than the 
spectrograph resolution, leading to typical errors of 25-50 m 
s"' per line. However, these calculations are specifically for 
the NIST SRM cells and custom-made longer pathlength CO 
cells (or more multiple passes through the same cell) can cre- 
ate deeper absorption lines if they are required. We speculate 
that significantly larger path lengths than provided by conven- 
tional cells may be achievable in the future by directly using 
long lengths of gas-fille d hoUow-core photonic crystal fibers 
now u nder development ( Benabid et al.l2005tlTuominen et al.) 
l2005h . These gas cells can also be used at higher (or lower) 
spectral resolutions than the R =50-70k we have considered 
here. To demonstrate this, we follow the prescription of 
Bouchy, Pepe & Queloz (2001) to calculate a Quality factor 
(Q) that is an estimate of the intrinsic radial velocity infor- 
mation contained in the spectrum. The Q factor depends on 
factors like the spectral richness of the wavelength region be- 
ing considered, the sharpness of the absorption lines, as well 
as the instrument profile. For a given wavelength region, the 
limiting radial velocity precision (ctrms) is given by 



where c is the velocity of light, Q the Quality factor, and 
the total number of photons collected. As expected, the radial 
velocity precision is proportional to ^/N. 

We calculate the Q factor for the HCN SRM for differ- 
ent spectral resolutions using the wavelength range 1524- 
1565nm. Figure [3] shows the Quality factor as a function of 
the spectral resolution. At low spectral resolutions all the in- 
dividual lines blur into each other, leading to very low ve- 
locity information content. At intermediate spectral resolu- 
tions (R=50-100k) the information content is a strong func- 
tion of spectral resolution, as the absorption lines, which have 
line widths of 7-12pm, are still unresolved. When the spec- 
tral resolution exceeds 200,000, the individual lines them- 
selves begin to be resolved and the Q factor begins to plateau, 
asymptotically approaching its finite value at very high res- 
olutions. Beyond a resolution of R=300k the lines are fully 
resolved, and the Q factor increase is minimal. Currently 
planned NIR high-resolution instruments all operate in the 
R=20-100k regime, and the HCN & C2H2 SRM cells are 
well-suited to these. The commercially available CO cells, 
however, have larger line-widths, making them less desir- 



able for the higher spectral resolutions. Custom cells, how- 
ever, can easily be manufactured. Other than high-resolution 
spectrographs, the absorption cells can be used in other NIR 
instruments. Planned H band spectrographs like the fiber- 
fed APOGEE instrument for SDSS-III (Allende-Prieto et al 
. 2008), or dispersed fixed-delay interferometer instruments 
(Guo et al. 2006) in the NIR, can use such cells for wave- 
length calibration. 

It is important to recognize that the absorption line spacings 
in these cells are sparse when compared to absorption line 
densities in typical M dwarf spectra. The high continuum S/N 
is absolutely necessary to ensure that the error in the wave- 
length calibration is smaller than the achievable precision on 
M dwarfs. This is one of the reasons we have discussed this 
approach only in the context of a fiber-fed instrument. If such 
cells were to be used for simultaneous calibration by passing 
starlight through them, then the achievable precision may be 
severely limited by the reference rather than the star Such 
an approach should be considered only when the need to cal- 
ibrate out instrument drifts is more important than achieving 
close to photon noise limited precision on the stellar target. 

5.3. Wavelength Accuracy & Stability: 

The vacuum wavelengths of the line centers for the species 
in the SRMs can, and have, been measured very accurately at 
low temperatures and pressures. For example, the measured 
vacuum wavelengths a t low pressure for the 3;/ lines of '^C'^0 
(Picqu e& Guelac hvili 1997) agree with the theoretical ly cal- 
culate d values from the HITRAN database (Rothm an et alJ 
l2005h to 0.02 pm (or < 4 m s"' at 1.6 um). For higher cell 
pressures, the lines begin to broaden and also shift due to in- 
teraction of molecules during elastic collisions. Accounting 
correctly for this pressure shift is the dominant source of un- 
certainty for determining the line centers of the absorption 
lines. This shift is accounted for by explicitly measuring it at 
various pressures and extrapolating a linear relationship. The 
resulting line centers for the higher pressure CO gas SRMS 
are certified by NIST to 0.4-0.7 pm (all 2a errors), HCN from 
0.04 - 0.24 pm an d C2H2 from 0.1-0.6 pm (Swann & Gilbei 
I2OOO. 2002. l2005h . So the line centers themselves are known 
to an accuracy of 4-60 m s"^ (la errors), quite comparable 
with the Th-Ar emission line errors quoted by PE83. A sig- 
nificant fraction of this error is actually the uncertainty in the 
pressure of the cells themselves, which may be possible to 
better constrain with custom-made cells. If a higher accuracy 
is really necessary, then high-resolution FTS spectra can be 
used to independently calibrate the cells. Instruments known 
to possess high stability over a few hours can also take the 
approach of Lovis & Pepe (2007) by acquiring many Th-Ar 
exposures to build up high S/N to determine a wavelength so- 
lution, and using that solution to determine the line centers of 
the absorption cells. 

Of perhaps more concern is the inherent stability of the ab- 
sorption lines themselves. A change in operating tempera- 
ture affects the pressure of the cell, which, due to the pressure 
shift, leads to a shift of the line center The wavelength shift 
AA due to pressure at temperatures T and T„i are related by 
(ISwann & GilbertlllOOOl) 

AMT) = AX{T,„)^{T/T,„). (2) 

A one degree change in temperature at a operating tempera- 
ture of 296 degrees leads to a 0.17% change in the pressure- 
induced wavelength shift. For the high pressure CO cells. 
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the pressure-induced shift is at most 3pm, and that one de- 
gree change leads to a line shift of ~ 0.005pm, which is a 
sub-ms"' effect. Given the fairly loose temperature require- 
ments the entire enclosure containing the cell assembly (Fig- 
ure |2| could be temperature stabilized to ±1C using heating 
elements running in a proportional-integral-derivative (PID) 
loop with temperature sensors. If necessary, temperature con- 
trol to ±0. 1 degrees is easily achieved if each cell is stabilized 
independently with a heating element wrapped around it, and 
can make temperature effects negligible. The absorption cells 
are relatively immune to any other environmental effect since 
they are contained in sealed glass tubes. Unlike Th-Ar lamps, 
absorption cells do not have a warm up period (typically 15-30 
minutes) where the lines are not stable enough for precision 
applications, or a limit on their operating lifetime. They are 
also passive, requiring only a light source for operation. 

5.4. Bright White Light Source 

The single-mode fiber-coupled gas cells makes for a com- 
pact and easily configurable set of wavelength references, but 
the narrow core and spatial filtering properties of such fibers 
make them inherently difficult to couple to traditional illumi- 
nation sources like quartz lamps. A number of white light 
sources coupled to single mode fibers are now commercially 
available for the telecom S, C and L bands (1400-1630nm), 
which are well matched to the spectral regime covered by 
these refereence cells. Alternatively, tunable diode lasers are 
also available that can be used to scan through the regions of 
interest at many times per second. Such sources may easily 
be adaptable for spectroscopic applications in astronomy. The 
science fiber coupled to the spectrograph needs to be multi- 
mode to couple effectively to the telescope. To ensure that the 
science and calibration fiber illuminate the spectrograph op- 
tics in the same way, the calibration fiber should be identical 
to the science fiber. Expanding the light from a single mode 
fiber to the correct focal ratio and coupling into the multi- 
mode science and calibration fibers poses no major challenge 
if this is deemed necessary. As mentioned before, the use of 
single-mode fibers is not mandatory if multi-pass configura- 
tions are not required. In such a case, a multi-mode fiber is 
adequate, still enabling the light to be passed through multi- 
ple cells. The penalty of this approach is the unwieldy 80 cm 
length of the CO cells, and the advantage is that a conven- 
tional quartz-lamp light source is quite sufficient. 

5.5. Isotopologues 

The NIST standards are deliberately designed to span as 
large a range as possible in the telecom bands. For such appli- 
cations, line densities are not as critical as for the broadband 
wavelength calibration application in astronomy. Available 
line densities and wavelength coverage can be increased by 
using isotopologues of these molecules, which exhibit sim- 
ilar levels of stability but have absorption lines at different 
wavelengths. For CO the NIST cells are specifically cho- 
sen to be i^ci^o and '^C'^O. The HITRAN database con- 
tains line transition parameters for 4 other CO isotopalogues: 
'^C'^O, '^ci^o, i^c'^O, and i^c'^^O. Figure H shows the ab- 
sorption lines for all 6 CO isotopologues listed above. To- 
gether, they substantially increase the line densities available 
for calibration in this region of the H band, and extend cov- 
erage to longer wavelengths. The filled circles in the figure 
correspond to zero-pressure line centers from the HITRAN 
database. Transmission values are all scaled to that of '^C'^O 
because the purpose of the figure is to illustrate the increase 



in density and coverage. Actual transmission is best derived 
from experimental data for such cells and the length required 
for gas cells of the four additional isotopologues may be dif- 
ferent than those of '^C'^O and '^C'^'O. Even higher line 
densities and coverage than shown here may be possible us- 
ing '"*CO, even though '"^C is unstable and undergoes jS de- 
cay to ''^N with a half-life of ~ 5730 years. Similarly the 
many known isotopalogues of '^C2H2 and HCN can also be 
used to generate additional absorption lines. While the use 
of isotopologues is attractive for increasing line densities, the 
pressure-shifted line centers for many of these species need 
to be determined before they can be used. Bootstrapping a 
wavelength solution off the known NIST calibrated lines is 
possible for very stable instruments, as described earlier. In 
general, however, the isotopologue gas cells should also be 
calibrated independently, like the NIST SRMS. 

5.6. Additional Molecular References in the NIR J & H 
bands 

Thus far we have primarily discussed the use of the NIST 
absorption gas cells and their isotopologues in the NIR H 
band. As discussed earlier, gas cells calibrations in this wave- 
length regime exist primarily due to the needs of the telecom- 
munication industry. Although such well-characterized cells 
are not commercially available in the J band, it is worth briefly 
discussing the molecules that are promising possibilities in 
this wavelength regime as well. Hydrogen Fluoride has a se- 
ries of sharp absorption lines in the 867-909nm and 1257- 
1340nm regions. Methane and water have lines in the J & H 
but both are also atmospheric species. D' Amato et al. (2008) 
have demonstrated that HCl, HBr and HI exhibit absorption 
lines in the 1.2-1.4/im region, though getting deep absorp- 
tion lines requires path lengths approaching 1 meter. Fiber 
fed gas cells similar to the NIST SRMs can be easily used to 
provide the path lengths necessary for moderately deep ab- 
sorption lines in such cases. HCL has a series of lines in the 
1 1 85 -1240nm region and 1 720- 1870nm regions, which well- 
complements the molecules we have discussed for the H band, 
and spans parts of the J band. 

5.7. Safety Considerations in Dealing with Gas Cells 

For practicality, the gas cells must not be toxic or lethal 
since a leak, or a breakage, would significantly impact oper- 
ations. Acetylene is not known to be toxic even if inhaled 
in large concentrations. Hydrogen Cyanide is lethal in large 
amounts since it inhibits enzymes in the electron transport 
chain of cells, thereby making normal cell functioning im- 
possible. The trace amounts of HCN found in the NIST 
SRM cells is quite safe. Even if the HCN contents on the 
entire SRM gas cell were to be respired- the increase in the 
cyanide content in blood would be minimal. Inhaled Carbon 
Monoxide is lethal in large concentrations since it binds to 
haemoglobin, preventing oxygen transport. However CO is 
also harmless in the small trace amounts found in the gas 
cells. Methane and water are commonly occurring atmo- 
spheric species and are not toxic. Hydrogen Fluoride is very 
corrosive and attacks glass. HF gas cells have to be specially 
designed with materials that do not react with the gas and the 
cells have to be handled carefully. Such HF cells have been in 
use for calibration for a number of years now. The gas cells 
we have explored here are all safe to use in a laboratory or a 
confined environment and breakages or leaks, if they happen, 
are not a cause for major concern. 
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TABLE 1 

Commercially Available NIST Wavelength Reference Gas Cells suitable for H band wavelength calibration. 



NIST SRM 


Species 


Range (nm) 


Pressure (Torr) 


Lines 


2517a 


C2H2 


1510-1540 


50 


56 


2519a 


h'^c'-'n 


1530-1560 


25 


51 


2514 


12C'60 


1560-1595 


1000 


41 


2515 


13C'60 


1595-1630 


1000 


41 


6. DISCUSSION 



We have explored the relative advantages of different wave- 
length calibration techniques for high-resolution fiber-fed 
NIR echelle spectrographs. We conclude that a series of com- 
mercially available absorption cell standards can be used to 
wavelength calibrate echelle data over a significant fraction of 
the H band, covering over 120nm with four gas cells. Some 
of these cells require long path lengths, but the use of single- 
mode fibers enables compact multi-pass configurations with 
small diameter cells that can easily be integrated into a cal- 
ibration unit. Although the absorption lines are very stable, 
their absolute line centers are not known a-priori to better 
than 4-60 m s~'. In principle, this uncertainty can be re- 
dressed by explicitly acquiring a high-resolution FTS spec- 
trum or measuring them with tunable diode lasers. In prac- 
tice, even this solution may not be necessary for most ap- 
plications. We hope to demonstrate the achievable precision 
of these cells with FIRST, a high-resolution silicon immer- 
sion grating bas ed instrument in the preliminary design stage 
jGe et al.ll2006l) . and to further explore their potential use as 



calibrators. We have also considered only the four commer- 
cially available NIST cells. By using additional gases, iso- 
topologues and cell lengths, one may be able to extend this 
technique to span larger regions of the H band with more and 
deeper absorption lines. Many NIR instruments in the de- 
sign or construction stage may benefit from using the calibra- 
tion approach we have outlined here. Such an approach may, 
in concert with Th-Ar lamp used simultaneously for regions 
without a absorption reference, be able to provide adequate 
wavelength calibration for most high precision applications 
until frequency-comb technology is mature, commonly avail- 
able, and less expensive. 
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NSF AST-0705139, NASA with grants NNX07AP14G and 
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Florida. 
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Fig. 1 . — Expected transmission for the NIST SRM cells when observed with a spectrograph with a spectral resolution of 50k. For HCN and C2H2, we have 
simulated the transmission at this resolution using available high resolution data of the NIST SRMs. The spectra shown are for a 5 cm long C2H2 cell and for a 
15 cm long HCN cell. For CO, the spectrum is simulated to resemble the spectrum observed by Swann & Gilbert (2002) using a O.OSnm bandwidth spectrum 
analyzer (effectively R=32k). The CO lines are expected to be a few percent deeper than shown here when observed at R=50k. While the length of the C^O cells 
is 20 cm, the light makes four passes through the cell for a total absorption path of 80 cm. 




Fig. 2. — A schematic showing the SRM gas cells, their lengths and an arrangement of cells that yields a transmission spectrum similar to that in Figure[T] 
The schematic illustrates the single-pass transmission for HCN and C2H2 gas cells as well as the 4-pass transmission in the 20 cm CO gas cell to yield a 80 cm 
absoiption path length. All cells are linked with optical fibers. For clarity we have shown the multi-pass arrangement as four separate beams, while in reality two 
sets of the beams are parallel. A full diagram of the CO SRM cell and holder can be found in Figure 2 of NIST publication 260-146 " 

''http://ts.mst.gov/MeasurementServices/ReferenceMateiials/upload/SP260-146.PDF 
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Fig. 3. — The Quality Factor (Q) for the HCN SRM as a function of the spectral resolution. The quality factor is an estimate of the radial velocity related 
information content in the spectra. The increase in the Q factor begins to plateau as the width of the resolution element becomes significantly smaller than the 
intrinsic hne-widths in the spectra. 
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Fig. 4. — Absorption spectra of isotopologues of CO that are present in the HITRAN database. Points coiTespond to the zero-pressure line centers listed in 
the database. Transmission spectra of all isotopologues has been scaled to that of ''c'^O for ease of comparison. The transmission values on the y-axis are for 
illustration only. 



